Hopping modulation in a one-dimensional Fermi-Hubbard Hamiltonian 
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We consider a strongly repulsive two-component Fermi gas in a one- dimensional (ID) optical 
lattice described in terms of a Hubbard Hamiltonian. We analyze the response of the system to a 
periodic modulation of the hopping amplitude in presence of large two body interaction. By (es- 
sentially) exact simulations of the time evolution, we find a non-trivial double occupancy frequency 
dependence. We show how the dependence relates to the spectral features of the system given by the 
Bethe ansatz. The discrete nature of the spectrum is clearly reflected in the double occupancy after 
long enough modulation time. We also discuss the implications of the ID results to experiments in 
higher dimensional systems. 



Ultracold atomic gases systems couple weakly with 
the surrounding environment and are highly controllable 
P) S S S]' therefore they offer excellent possibihties 
to investigate the dynamics of strongly correlated quan- 
tum many-body systems. Much attention has been re- 
cently devoted to the study of the dynamical proper- 
ties both from the experimental 0, [6[ and theoretical 
0,[l,[l,[lOl, JJa] point of view. Especially, one-dimensional 
(ID) systems, accessible by experiments and theoretically 
exactly solvable in some cases, can be used to obtain thor- 
ough understanding of the many-body ground state and 
the dynamics. In this letter, we present an (essentially) 
exact time-evolving block decimation TEBD simulations 
of the dynamics of a repulsively interacting ID system 
and reveal a non-trivial time-dependence which we ex- 
plain using the Bethe ansatz (BA). We extend the anal- 
ysis also to the harmonically trapped case essential for 
ultracold gases experiments. In higher dimensions, the 
relevance of Mott and antiferromagnet (AF) physics in 
connection to high- Tc superconductivity (see e.g. [l3]), 
suggests that the investigation of the equivalent systems 
in the framework of ultracold gases, especially in two di- 
mensions, may shed new light on high- Tc superconductor 
physics. Our results are relevant for such experiments by 
showing — with analysis that does not assume mean-field 
approximation nor linear response — how the discrete 
nature of the spectrum is reflected in the dynamics of a 
lattice modulation experiment. 

We examine the dynamical properties of a two-species 
ultracold atomic gas loaded in a ID optical lattice both 
for open boundary conditions (obc) and in presence of 
parabolic confinement. In particular we perform an (es- 
sentially) exact numerical simulation of this system when 
a periodic lattice modulation is applied. To this end we 
consider the ID Hubbard Hamiltonian, in presence of an 
external parabolic confining potential 
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h.c. Hi, 



(1) 



U / J ^ 1, the only effect of the lattice modulation is 
a modulation of the ho ppi ng amplitude J, the effect on 
U being negligible (see [ll|). 

We focus on the double occupancy (d.o.) expectation 
value < n^^^nil > , when a (small) periodic modulation 
5Jsa\(ujt) of the hopping amplitude J is applied for a 
given time to the ground state of the Hamiltonian given 
in Eq. ^ for two different situations: J7/ J = 60, L = 20, 
particle number Np = 12 and U/J = 20, L = 40, Np = 24 
respectivel y. O ur simulation is performed with a TEBD 
algorithm [l3|, [ij] , both for the ground-state calculation 
(imaginary-time evolution) and the real-time evolution. 
The numerical results show that the ground state is con- 
stituted by a central Mott region with one atom per site, 
surrounded by two small metallic (Luttinger liquid) re- 
gions where the filling is less than one. In order to avoid 
finite-size effects, we have considered a lattice size ex- 
ceeding the actual extent of the atomic cloud by a few 
lattice sites. Heuristically, the Luttinger liquid phase cor- 
responds to the regions where < nf > — < rii const, 
see [l5|. Moreover, the static structure factor S{q) in the 



where Hj = - cL^i+i cr ^ "--t^-, ^J-int 

U^f nn^riii, Vi = — io)^ , J is the hopping am- 
plitude, and U the on-site interaction. In the limit 



ground state of the finite systems here considered ex- 
hibits the same qualitative features of the ID Heisenberg 
AF chain (i.e. a slow decaying peak centered around 
q — tt), as expected in the U ^ J limit of the Hubbard 
Hamiltonian [16|. 

For an infinite chain at half fiUing, if hopping and 
parabolic confinement are suppressed, the (highly degen- 
erate) first excited state is represented by a site with 
an empty site and a doubly occupied one (henceforth 
particle- hole excitation). The energy gap between this 
state and the ground state is equal to U. To investigate 
the spectral properties of the system when [/ ^ J, it 
seems then natural to choose w ~ [/ in SJsin{Lot), with 
5J / J = 0.1 throughout the paper. 

In Fig. [1] (inset) we show the d.o. as a function 
of frequency for short times, from which one notices 
that a broad peak appears as in correspondence of the 
value = [/, consistently with the value of the gap of 
the particle-hole excitation. However, for larger times a 
richer structure appears (Fig. [T]). From Fig. [21 it is clear 
how it is possible to distinguish between short- [t < t*) 
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FIG. 1: Double occupancy < ni^^ rii^i > as a function of 
frequency for long {t = 10 main panel) and short {t = 0.5 
inset) times. A single broad peak appears for LojU ~ 1 (L = 
20, C//J = 60, n/J = 0.1). 
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FIG. 2: (color online) Double occupancy J]]^ < ni.in^^i > 
as a function of time, for two different frequency values cor- 
responding to a minimum {uo/U — 0.978) and a maximum 
{i^/U = 0.988) of the plot in Fig. [T](same parameters as Fig. 



and long- (t > t*) time behavior, the exphcit value of the 
threshold t* will be derived below. 

The key idea in understanding the time dependence 
of < niijUii^ > is that the resolution of degeneracy de- 
pends on the modulation timescale. For t < t* it is 
possible to consider the particle-hole excited states as 
quasi-degenerate, hence the transition probability be- 
tween the ground state and the quasi-continuum (cen- 
tered around U) of excited states is given by P{t) oc 
5.P sin^ [{U - hiu) t/ 2h] {U - hwf. The threshold time 
is then calculated as the time when the internal struc- 
ture of the excited states band becomes visible, namely 
t* = 2tt/ where /S.E is the energy difference in the 
quasi continuum of excited states which, in our case, will 
be determined by BA. For t > t*, the transition proba- 
bility is given by 



IK 



(2) 

where En represents the energy of one of the quasi- 
degenerate excited states. The numerical results plot- 
ted in Fig. [5] show how the transition between these 
two regimes is quite abrupt, allowing an easy compari- 
son with the theoretical value of t*. Eq. ^ explains 
the piecewise linear time dependence of d.o., which is 
induced by the different effective density of states p{En) 
and perturbation matrix elements |Ki,gs(<5J)P before and 
after t* . In particular, in order to have a peak in the 
d.o. spectrum, both the perturbation matrix element 
and the density of states at the appropriate energy must 



be different from zero, suggesting that the d.o. spec- 
trum exhibits selection rules in its long-time peak dis- 
tribution. In the 2-site problem, it is possible to see 
explicitly how the ground state with two particles, rep- 
resented by singlet state 15 >= (| Tji> — I ijT>)i is 
coupled by the hopping modulation to the state |Z)+ >= 
75 (I Ti, > +|0, ti>) only, and not to the state \D^ >= 

^(1 n,o>-|o,n>). 

The ID nature of the problem allows some insight on 
the peak position deriving from the exact (BA) solution 
of the ID Hubbard Hamiltonian. As a first step we make 
contact between BA and the numerical solution of a small 
linear chain with open boundary conditions {Vi = 0) and 
U/J — > oo{U/J = 500). To this end we consider the 
BA equations for an open Hubbard chain (see [l3| and 
EPAPS supplementary material [l^]). The solution of 
BA equations with respect to the charge momenta kj 
and spin rapidities Ac allows to determine energy and 
momentum eigenstates whose values can be expressed in 
terms of ki as 
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Along the lines of the derivation by Ogata and Shiba [Iq 
applied to the case of OBC, Bethe-ansatz equations lead 
to the simple relation kj — n Ij / (L + I), in the limit 
U/J — !■ oo. The latter expression is particularly relevant 
since it describes the spectrum of the considered Hamil- 
tonian in terms of spinless fermionic particles. This re- 
sult can be interpreted within the general framework of 
spin-charge separation of excitations in ID systems (see 
e.g. In particular the ground state for iV = L 

is obtained when Ij — 1 ... TV. Its energy is given by 
E = — 2J^^^-^ cos{kj) and the Fermi quasi-momentum 
is given by hp = n L/{L + 1). The excitations of the 
system can be described in terms of particle-hole excita- 
tions of this system. Hence, since the first Hubbard band 
is full, 

AS ^ AS + [/ = -2J(cos(A:p) - cos(fc/,)) + [/, (4) 



where Ep — 



'2 J cos kp + U corresponds to the energy 



of an extra particle added in the second Hubbard band, 
and 2J cos{kh) the energy of the hole in the first Hub- 
bard band. Since having a particle and a hole in the 
same momentum state would not contribute to the in- 
crease of the interaction energy by U, these states must 
not be considered in the calculation of the energy level 
structure around U. We also note that in the thermody- 
namic limit L — ^ oo, the discrete energy-level structure 
becomes a continuous band of width 8J/U, as it is possi- 
ble to deduce from Eq.jjl]). Our simulations are at T = 0, 
but the features in the spectra are expected to be visible 
if the system is in the Mott state and T < AE. 
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FIG. 3: (color online) d.o. spectrum for U/J = 500 and 
L = 4,6, 8, 10 (top- left to bottom right). Vertical lines repre- 
sent the value AE for various kp and kh- Red continuous lines 
represent the values listed in Table |I] Blue dashed lines rep- 
resent transitions between states of even parity (i.e. ip and ih 
even), which do not seem to correspond to a peak in the d.o. 
spectrum. Grey dotted lines represent states with ip = 2n + l 
or ih = 2n + 1, with n € Z . Green dash dotted lines rep- 
resent a particular condition of coincidence of the frequency 
between a "red line" transition and a "gray line" transition 
with ip — 2n + l, ih = Iri + 1 and ip — ih ~ 2n" (see Table H}. 
For ease of reading, the right part of the spectrum has been 
omitted due to its symmetry around u} = 1 (c.f. the slight 
asymmetry in the trapped case, Fig. [T]). x axis: uj/U, y axis: 
< nifmi >). 



TABLE I: Theory /numerics comparison 
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We have made a connection to the particle-hole ex- 
citation spectrum by performing hopping modulation on 
systems with open boundary conditions (i.e. no parabolic 
confinement), with U/J — 500 and L — A, 6, 8, 10. As 
seen in Fig. [3l the d.o. peaks correspond to specific 
particle- hole excitations, supporting the explanation of 
the d.o. spectrum in terms of selection rules (see TableHl. 
Even if the simulations in presence of a parabolic con- 



fining potential with U / J = 20, 60 cannot be described 
exactly within the above formalism, we suggest that it is 
possible to make contact between a description in terms 
of spinless fermions and the numerical results, the only 
effect of parabolic confinement being the different spec- 
trum in the Mott phase. This hypothesis is justified by 
the nature of the metallic phase in the limit U/ J ^ 1 
(gapless spectrum, small spatial extent) and confirmed by 
numerical evidence. As in the case where the parabolic 
confining potential is absent, we assume that the spec- 
trum of the system can be described in terms of particle- 
hole excitations of spinless fermions moving in a lattice, 
in presence of a global parabolic confinement. 

In [2^ , it has been shown that, for 4 J > Q, the approx- 
imate description of the single-particle spectrum can be 
carried out in terms of low-energy excitations with quan- 
tum number n < ric, Uc = 2||-\/2jT2|| — 1, where ||a;|| is the 
integer closest to and high-energy excitations with 
n > Tic. li n < He, the single particle energy is given by 



Eo = 2\/jT7(n 



32 



(2n+ 1)2 + 1 



1/2)- 
(2n+ 1)3 + 3(2n+ 1) 



(5) 



where Eq is a constant energy term. If, on the other hand, 
n > Tic, for the high energy modes the energy eigenvalues 
are given by 



E„ 



En 



2J 



(2z)2 



r 



(6) 



representing states close to position eigenstates. 

In the trapped system, performing a calculation anal- 
ogous to the case with open boundary conditions, it is 
now possible to compare the particle-hole excitation spec- 
trum to the d.o. spectrum obtained numerically (see Fig. 
2]). It must be noted that, in Fig. 31 we have shifted 
{Auj/U ~ 0.667) the energy spectrum calculated in Eqs. 
([5]), (HI), in order to obtain the superposition with the 
d.o. spectrum. The shift is due to finite value of U. This 
effect can be explicitly calculated in the two-site case, 
where the energy-level splitting can be obtained analyt- 
ically. In the case U/J = 60, the d.o. spectrum shows 
good agreement with the spectrum from Eqs. (O and (O, 
when we enforce the condition Up = even, = even. 

Our analysis of the dynamical properties of a ID Hub- 
bard chain is in close relation to the experimental [5,] and 
theoretical ll| investigations of Mott and AF phases in 
optical lattices. Even if in our simulations there is no 
proper AF phase, due to the dimensionality of the sys- 
tem, it is possible to relate the results obtained here con- 
cernin g th e peaks in the d.o. spectrum to those appear- 
ing in [ll|. In presence of AF ordering, an energy gap 
cx J^/U develops. In our case, the separation between 
the different energy levels, which mimics the effects of 
an energy gap, is induced either by the parabolic con- 
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FIG. 4: (color online) Trapped system: energy levels (contin- 
uous red lines correspond to even values of ip, ih) for L — 20, 
Ntot = 12, U/J = 60, n/J = 0.1. The peaks at the edge 
of the diagram do not correspond to any peak due to the 
approximate description of the energy levels. 



lation, we have analyzed the properties of a fermionic gas 
in a ID optical lattice, in presence of parabolic confine- 
ment. The numerical results show that, for a sufficiently 
long time, a nontrivial peak structure appears in the d.o. 
spectrum. The peak structure has been qualitatively jus- 
tified in terms of the BA solution of a chain with open 
boundary conditions, in the limit of U/J — + oo. The re- 
sults we have obtained, while not being a direct evidence 
of what should be observed in a modulation experiment 
in 3D [1| when T < TnooI 21 \ , due to the different nature 



of the gaps in 3D and our ID systems, explain in detail 
how the discrete nature of the spectrum is reflected in 
the lattice modulation experiment. The results thus sug- 
gest that the gap in the AF phase can be observed via a 
nontrivial peak structure. 

This work was supported by Academy of Finland 
and EuroQUAM/FerMix (Projects No. 213362, 217041, 
217043, and 210953) and conducted as part of a EURYI 
scheme grant |23j]. 



FIG. 5: (color online) Trapped system: energy levels (same 
definitions as in Fig. Ell L = 40, iV = 24, f/ = 20, Q/J = 0.1 
(right), when values of ih and ip up to 20 are considered. 



finement or by finite size effects. The qualitative sim- 



ilarity between the spectrum obtained here and in 11 | 
confirms the relationship between the presence of finite 
energy gaps and the peaks in the d.o. spectrum. In 
our case, especially for the obc case and for U = 60, we 
are able to suggest how the excitation spectrum together 
with the nature of the lattice modulation determines the 
structure of the peaks. While [ll[ considers a higher di- 
mensional system thus closer to the present experiments, 
our study is complementary in the sense that it does not 
assume mean-field approximation or linear response. 

In addition, when performing a d.o. modulation exper- 
iment aiming at the observation of the spectrum depicted 
in Fig. [U the effect of a finite modulation time must be 
taken into account. If we consider experimentally rel- 
evant values (see [5]) of U, J and the modulation time 
tmod {U ~ 5kHz, J ~ O.lkHz, t^od = 50/C/ = 0.01 s ), we 
obtain t* ~ 2nU / J'^ ~ 3.14 s. Since this is experimen- 
tally rather long time, the estimate suggests that smaller 
values of U / ,P might be necessary to observe the AF gap. 
In this calculation of t* we have considered /S.af — J'^ /U . 

In summary, with the aid of a TEBD numerical simu- 



[1 
[2: 

[s: 
[4: 
[s: 
[6: 

[r 

[9: 

[10 

[11 

[12 

[13 
[14 
[15 
[16 

[17 
[18 
[19 

[20 
[21 
[22; 



Electronic address: paivi.torma@hut.fi 
D. Jaksch and P. ZoUer, Ann. of Phys. 315, 52 (2005). 
Chin et al., arXiv:0812.1496 (2008). 
Jaksch et al., Phys. Rev. Lett. 81, 3108 (1998). 
Peil et al., Phys. Rev. A 67, 051603(R) (2003). 
Jordens et al.. Nature 455, 204 (2008). 
Schneider et al.. Science 322, 1520 (2008). 
KoUath et al., Phys. Rev. A 74, 041604(R) (2006). 
Winkler et al.. Nature 441, 853 (2006). 
Cramer et al., Phys. Rev. Lett. 100, 030602 (2008). 
Flesch et al., Phys. Rev. A 78, 033608 (2008). 
Sensarma et al., 0902.2586 (2009). 
L. Hur and T. M. Rice, Ann. Phys. (N.Y.) 324, 1452 
(2009). 

C. Vidal, Phys. Rev. Lett. 91, 147902 (2003). 

G. Vidal, Phys. Rev. Lett. 93, 040502 (2004). 

M. Rigol et al., Phys. Rev. Lett. 91, 130403 (2003). 
F. H. L. Essler et al.. The One-Dimensional Hubbard 
Model (Cambridge University Press, 2005). 

H. Asakawa and M. Suzuki, JPA 29, 225 (1996). 

M. Ogata and H. Shiba, Phys. Rev. B 41, 2326 (1990). 
T. Giamarchi, Quantum Physics in One Dimension (Ox- 
ford University Press, 2003). 

A. M. Rey et al., Phys. Rev. A 72, 033616 (2005). 
J. Neel, Ann. Phys. (Paris) 5, 232 (1936). 
See EPAPS Document No. E-PRLTAO- 103-077934 for 
supplementary material on Bethe-ansatz equations for 
open-boundary conditions. For more information on 



EPAPS, see http://www.aip.org/pubservs/epaps.html 



[23] www.esf.org/euryi. 



